Dissipation effects in the ratchetlike Fermi acceleration 
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Ac driven asymmetric pulses can be used to control the Fermi acceleration between three different 
motions, A : the accelerated mode; D : the decelerated mode and H : the hyperaccelerated mode. In 
this work we show that dissipation strongly affects the particles velocity, reducing the possibility 
for an accurate control of the dynamics. The saturation time, where the mean velocity starts to be 
constant due to dissipation, decays with a power law ~ 7 _/3 , where 7 is the dissipation parameter 
and j3 is close to 1. The value of the saturated mean velocity also decays with a power law with 
exponent j3 ~ 0.6 for the case H, and j3 ~ 0.3 for the case A. In the case D this velocity is almost 
constant for small dissipations. 
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Fermi acceleration is a to pic which got attention in various areas of physics, ran ging from nonlinear physics [H, 0, S 
C*~) ! 0, H, 0, Ilk atom optics [il. I id. [ill], plasma physics 12, 3] to astrophysics (14 Fl5l.Tl6lf After the first model proposed 



by Fermi [171 ]. essentially two different versions became common in the literature. In the first one, the Fermi-Ulam 
(FU) model, a bouncing particle moves between a fixed surface and a parallel oscillating surface [3]. In this case the 
regular islands in the phase space prevent the Fermi acceleration. A simplified version of this model was proposed 
\ to improve simulations [lj, called the static wall model. It essentially ignores the displacement of the moving wall 
but keeps the information for the momentum transfer as the wall was oscillating. The dynamics of the static model 
was studied in different aspects [H, 0, 0, 0] and for different models [f| , and the relevant result for the purpose of the 
present work is that invariant curves in the phase space, found for higher velocities, prevent the particle to increase its 
kinetic energy without bounds. Recently the hopping wall approximation was proposed 0,0] which takes into account 
the effect of the wall displacement, and allows the analytical estimation of the particle mean velocity. Compared to 
^ I the simplified static model, the particle acceleration is enhanced. The second kind of Fermi accelerated model was 
proposed in 1977 by Pustyl'nikov (l9j |. who considered a particle on a periodically oscillating horizontal surface in the 
\^ \ presence of a gravitational field. Different from the FU model, for some initial conditions and control parameters the 

■ particle energy can grow indefinitely. 

j. I I In this paper we analyze the effect of dissipation in the simplified FU model when an ac driven asymmetric pulse 

■ controls the Fermi acceleration (deceleration) Small dissipation is inevitable in real systems and its influence 



Q\ \ on the dynamics of conservative systems is of great interest [2ll . [22|, [23| since elliptic periodic orbits become small 
sinks and attractors start to exists 0, [25| . In the context of Fermi models dissipation effects have been inserted in 
two ways: frictional force [2(| and inelastic collisions at the walls 27, 28|, 29[. Here we use the second approach and 



" ^ ■ the pulse is a deformed sawtooth driving law for the moving wall. This Ratchetlike pulse differs from the ac driven 
asymmetric pulses (symmetric sawtooth) used for the Fermi acceleration in the early work of Lichtenberg et. al [lj 
and proposed recently to control the motion of magnetic flux quanta [3(| and to analyse the relative efficiency of 
mechanism leading to increased acceleration in the hopping wall approximation Q . In the simplified Fermi model [H 
the particle is free to move between the elastic impacts with the walls. Consider that the moving wall oscillates 
between two extrema with amplitude vo. The gravitational force is considered zero. The system is described by a 
two-dimensional map M 1 ^ 2 )(V n , </>„) = (V^+i, <p n +i) which gives, respectively, the velocity of the particle, and the 
phase of the moving wall, immediately after the particle suffers a collision with the wall. Considering dimensionless 
variables the dissipative FU map with the deformed sawtooth wall is written as 
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for (j> n > 7/1, where n is the iteration number and /j is the maximum distance between the walls. For 7 ^ wc 
have dissipation effects. Since in this simplified model the displacement of the moving wall is ignored, the modulus 
function is used to avoid errors due to successive collisions which may occur in the original model. In other words, 
if after a collision with the wall the particle continues to have a negative velocity (a successive collision will occur in 
the original model), the particle moves beyond the wall. The modulus for the velocity injects the particle back and 
fixes the problem. 

The time asymmetry of the oscillating wall in Eqs. (TJ]) and is controlled by varying the parameters (771, 772). The 
deformed sawtooth (Ratchetlike) is obtained when 771 ^ 772. Figures [Ha)- (c) show the time behaviour of the oscillating 
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FIG. 1: (a)-(c) The shape of the pulses used in the simulations. The deformed sawtooth effect is obtained when An = 772—771 7^ 
0.0. For (a) An = we have the symmetric accelerated case A, (b) A77 < we have the deformed decelerated case D and (c) 
An > the deformed hyperaccelerated case H. 

wall (the pulse) for different values of the asymmetry A77 = 772 — 771. The deformed sawtooth pulse is obtained when 
A77 ^ 0.0. Such pulses can be easily obtained from pulse generators. For the case of no dissipation three different 
motions were obtained and controlled (2p| . They are A : for A77 = (symmetric sawtooth) the accelerated case [see 
Fig. [UJa)]; D : for A77 < 0(asymmetric sawtooth) the decelerated case [Fig. [5fa)l; and H : for A77 > (asymmetric 
sawtooth) the hyperaccelerated case [Fig. B(a)]. This classification, different from [3l|, is based on how fast the average 
velocity grows or decreases. At next we discuss separately the effects of dissipation in each case. To do this we show 
the corresponding phase space dynamics and determine the mean particles velocity at a given time n from 



t=o j=l 

where the index i refers to the ith iteration of the sample j, and £ is the number of initial conditions. We iterate the 
map ([T]) or for times n = 1 x I0 8 and 3000 initial conditions in the interval < (f> < i]i + 772 and < V < 10 -3 . 
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FIG. 2: Evolution of 500 chaotic orbits on the phase space V x <j> for the parameters n — 10, vo = 0.2 for case A (An = 
0.0, r\\ — r\2 — 1.00). For Fig. (a) 7 = 0.0 and (b) 7 = 10~ 6 . The mean values of the velocity, calculated over 3000 trajectories, 
are shown in (c)-(f) and are related to the cases 7 = 0.0, 10~ 6 , 10~ 5 , 10~ 4 , 10 _1 . 



First we consider the accelerated mode which is similar to the simplified model studied by Lichtcnberg and Lieber- 
man [lj,|2|, where the harmonic force was considered. This case is shown in Fig. [5J As the particle velocity increases, 
regular island are observed and (V) increases slowly until ~ 10 [Fig. Ufa) and (c)]. The regular islands prevent the 
particle velocity to increase very fast. We mention that all initial conditions start inside the chaotic region at low 
velocities. The growth rate of (V) depends on the number of regular islands inside the phase space. For very small 
dissipations (7 = 10 -6 ) the regular islands are transformed into sinks, as can be observed by comparing Fig. Ela) 
and Fig. [U(b) (observe the dark regions), and attract the chaotic trajectories which pass nearby. As a consequence 
the mean velocity cannot increase as before (compare Fig. [5](c) and (d)). This effect increases for higher values of 
the dissipation parameter, as can be observed in Figs. [DJe)-(f). Therefore in case A the dissipation decreases the 
mean velocity (to a constant plateau) for any values of 7 when compared to the non-dissipative case. There are two 
additional information we can get from Figs. (2^d)-(f), the saturation time (t s ), where the mean velocity starts to be 
constant in time and the saturated mean velocity (v s ). Plotting both quantities as a function of 7 in the interval 
10 -6 < 7 < 10 _1 we obtain the decay t c ~ 7 -10 and v c ~ 7 -0 3 . 
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FIG. 3: Evolution of 500 chaotic orbits on the phase space V x <j> for the parameters /j, — 10, vo = 0.2 for case H (A77 = 
0.01,7)1 = 1.00,772 = 1.01). For Fig. (a) 7 = 0.0 and (b) 7 = 10" 6 . The mean values of the velocity, calculated over 3000 
trajectories are shown in (c)-(f) and are related to the cases 7 = 0.0, 10~ e , 10 -5 , 10~ 4 , 10 _1 . 

The second case to discuss is the dissipation in the hyperaccelerated case H [see Fig. EKa)]. For a very small 
asymmetry An = 0.01 the phase space is totally filled and no regular islands are observed [Fig.|3fa)]. As shown in [20l |. 
the Lyapunov exponents approach zero for this case. The corresponding (V) increases very fast until (V) ~ 3 x 10 4 [see 
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Fig- Etc)], showing that the accelerated mode is enhanced when compared to the case A. By adding small dissipation 
(7 = 10~ 6 ) into the system we observe in Fig. GUb) that again many sinks appear. In this case however, since no 
regular islands were observed for 7 = 0, we are not able to relate the sinks to regular islands. On the other hand, 
we can say that the dissipation is able to find attracting points which were not visible in the conservative limit. Also 
here small dissipation prevents the hyperacceleration of particles. The mean particles velocity also decreases (to a 
constant value) due to the small sinks which appear in the phase space shown in Fig. E^b). For higher values of the 
dissipation parameter the maximal mean velocity decreases. Observe Fig. Etd)-(f) for 7 = 0.0, 10~ 6 , 10~ 5 , 10~ 4 , 10 _1 . 
For the hyperaccelerated case we found that the saturated quantities decay like t c ~ 7" 11 and v c ~ 7~ 6 . 

Now we discuss the decelerated case D (A?; = —0.01). Here many regular islands appear in phase space [see 
Fig-Et a )] which again prevent the acceleration to increase without bounds, as in case A. The regular islands arc very 
large so that the maximal instantaneous velocity is around ~ 2.5, while for case A it was ~ 70.0. The corresponding 
(V) remains here almost constant [see Fig. HJc)]. As in the other cases, when small dissipation is added (7 = 10 -6 ) 
into the system, regular islands are transformed into sinks. However here the mean velocity increases [see FigfJJd)]. 
This increasing in the mean velocity is easy to explain. Since the first (from below) four islands from the conservative 
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FIG. 4: Evolution of 500 chaotic orbits on the phase space V x cf) for the parameters fi = 10, vq — 0.2 for case D (An = 
-0.01, 771 = 1.01, rj2 = 1.00). For Fig. (a) 7 = 0.0 and (b) 7 = 10~ 6 . The mean values of the velocity, calculated over 3000 
trajectories, are shown in (c)-(f) and are related to the cases 7 = 0.0, 10" 6 , 10" 5 , 10~ 4 , 10~ 3 , 10~ 2 , 10 _1 . 

system are very large (see FigfJJa)), almost all initial conditions cannot reach velocities higher than <~ 1.1. It is a kind 
of upper partial barrier for the velocities. However, when small dissipation is added, these four island are transformed 
into sinks so that the chaotic trajectories can penetrate the islands and the upper barrier does not exist anymore. 
This makes the mean velocity to increase because the position of the sinks are higher (in velocity) than the upper 
barrier. Such increase in the mean velocity when dissipation is present was also observed in (25l |. For the decelerated 
case we found for the saturated time t c ~ 7 -1 ' 2 while v c is constant in the interval 10 -6 < 7 < 10~ 2 and decreases 
for 7 = 0.1. 

To conclude, for a long time the Fermi acceleration has been studied in different models and applications [l], @, 0, 0, 
IE 0, S, IE EE ED, El III EI EE EE E3, EE EE EE , but just recently [U a deformed sawtooth (Ratchetlike) pulse 
was proposed to control the Fermi de(acceleration). Changing the asymmetry parameter from the Ratchctlikc pulse 
it was possible to get Fermi hyperacceleration and deceleration. By switching the pulse between hyperacceleration 
and deceleration modes an accurate control of the particles velocity was achieved. Beside the remarkable control of 
velocities obtained in the dissipation free problem, we observe here that dissipation effects transform regular islands 
into sinks, change the mean velocity and reduce the ability to control accurately the particle dynamics. The saturation 
time decays with a power law t c ~ 7 , where (3 is close to 1.0 (±0.01) for all considered cases. On the other hand, 
the saturated mean velocity also obeys the power law decay, but with exponents (3 ~ 0.6 for the case H, and (3 ~ 0.3 
for the case A. In the decelerated case the saturated velocity is almost constant for small dissipations. Our results 
motivate further analysis related to the control of Fermi de(acceleration) for larger Ratchetlike asymmetries of the ac 
driven pulse in the presence of dissipation [32| and in the hopping wall approximation. In the later case we expect to 
obtain a better efficiency to control the particles velocity since the the energy gain per collision is not underestimated 
as in the present model. In addition, it would be interesting to analyse the implementation of the Ratchctlikc ac 
driven pulse applied to the walls of quasi one-dimensional billiards which are coupled to thermal baths at different 
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temperatures. In such cases it is desirable to optimize and achieve the control of directed heat conduction [331 ] . 
which could leave to an important increase of the thermal efficiency in physical devices such as rectifiers and thermal 



transistors 



34] 
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